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Abstract 

We extend the classical approach of the ^-separation of the Laplace 
equation = (as a null eigenvalue problem) to the general steady 
state Schrodinger equation including cases with a scalar potential V and 
the energy is a fixed constant. 
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1 Introduction 

Let us consider the Schrodinger equation 

_| A ^ + (y-^ = 0, (1.1) 

on a Riemannian manifold Q with metric tensor G of any signature. According 
to Moon and Spencer [T5J HH1 UZj , if the ansatz 

</> = RlliM<t), 

where R is a suitable function of the coordinates and each fa is a function of the 
corresponding coordinate q 1 only, permits to split equation ()1.1|) into n separated 
ODE's, then this equation is said to be R-separable. We are interested in the 
following question: under which conditions on the metric G, the potential V and 
the value of the energy E is this ansatz successful? 

Usually, in the literature [TSJ EH H3 HUJ HSl El U\, authors require that in- 
separable solutions for Schrodinger or Heimhoitz equations exist for all possible 
eigenvalues E. In the present paper we study the more general theory in which 
-R-separation occurs for a single eigenvalue E, as in the Laplace equation, where 
the only eigenvalue is E = 0. We call this theory fixed energy R-separation (FER- 
separation). As a result, we shall see that FER-separation enlarges the number 
of classes of separable coordinates. First of all we analyse the definition of FER- 
separation, showing that the functions fa must depend on a suitable number of 
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parameters (c a ) in order to satisfy a completeness condition (Definitional- Ac- 
cording to this definition, we will provide a criterion for testing if FER-separation 
occurs in a given coordinate system (q l ) (Theorem l30|) : we shall see that a nec- 
essary condition is that coordinates be orthogonal and conformal separable (as 
defined in Then, we give necessary intrinsic conditions for the existence of 
FER-separable coordinates based on conformal Killing tensors (Proposition EZj)- 
Furthermore, ODEs whose integrations provide /^-separated solutions of are 
deduced and the associated second order commuting operators are determined. 
In such a way, we see that our definition of /^-separation is equivalent, when 
restricted to the Laplace equation, to the definition of orthogonal separable co- 
ordinates given by Kalnins and Miller 9 involving a set of operators constructed 
from a Stackel matrix. Indeed, the R-separation of Laplace equation is charac- 
terized in [5| by means of conformal symmetry operators. We remark that in our 
approach, both the orthogonality and the existence of a related Stackel matrix 
are proved and not assumed. 

Also the intermediate results, the "tools" used for proving Theorem |3U1 ap- 
pear to be quite interesting by their own: they are (i) a geometric interpretation 
of the additive separation of variables in a single null PDE of any order, (ii) 
three criteria for additive separation of such a PDE, which extend those given 
for the null Hamilton- Jacobi equation in jlj, both exposed in Sectional This 
geometrical approach is quite general and easily adapted for dealing with fixed 
energy separation and FER-separation of other PDE's of physical interest, like 
time-dependent Schrodinger equation or Schrodinger equation with a vector po- 
tential. For instance, as a straightforward application of the method, we get 
that, on an Einstein manifold, ordinary separation for a single value of E of 
only occurs in the coordinates allowing separation for all admissible values 
of E (this is well-known for the Laplace equation, which on Einstein manifolds 
is separable in the same coordinates separating Helmholtz equation). In other 
words, studying the fixed energy case of ordinary separation is not so interesting, 
since in the most common cases we do not get any new coordinate system, while 
for FER-separation the class of possible coordinates is definitively bigger. This 
is explain the importance of FER-separation, whose theory, even if applied to 
Laplace equation till the 19-th century, needs a precise mathematical foundation. 

2 From multiplicative R-separation to additive separa- 



The link between .R-separation of (jl.lj) and additive separation of a suitable PDE 
depending on G, V, R, and E is obtained via the usual (see for instance [TU] ) 
substitutions, 



tion 



u — h\(j) 
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ip = R4>. (2.2) 

We recall that the Laplace-Beltrami operator A = g^VjVj on a n-dimensional 
Riemannian or pseudo-Riemannian manifold (Q, G) with contravariant metric 
tensor G = in local coordinates has the form 

Aij = flr^'SgV - T h d h ij. (2.3) 

where T h are the contracted Christoffel symbols, defined as 

T h = g^ p T k = g kh T h . (2.4) 

Proposition 1. There is a one to one correspondence between the solutions of 
(ji.ij) of the form ip = R\\ i ^{q 1 ) and the additively separated solutions u = ln0 
of 

gV UiUj + g tt u u - fX + §E-U = 0, (2.5) 
where Ui = diu, Uu = d\u, U is the modified potential 

and 

P = g^(T j -2d j \nR). (2.7) 

Proof. By inserting (j2.2j) in (jl.lj) for a real function R ^ on Q, we see that ()2.2I 
is a solution of (jl.lj) if and only if satisfies 

RA(j) + 2Vi? • V0 + (%R(E -V) + AR) = 0. (2.8) 

which is equivalent to 

A0 + 2Vlni?- V0 + (Jr£- £7)0 = 0, (2.9) 

where C7 is given by (j2.6j) . Then, let us perform the substitution (j2.1j) and let 

us denote the partial derivatives of u by = diU, = dfjU, Uijk = dfj k u, etc. 
Since u is additively separated, we get = 6 fan. Moreover, by (j2.1j) we have 

<%0 = «i0, dlcj) = (uu + m-)0, dijcp = UiUjcj) (i ^ j). (2.10) 

Therefore, by inserting (j2~lT)j) in (j2~9"j) . due to (|2~lj) . we get 

{g ij u iUj + g u u u - fX + Jr£ - C/)0 = 0. 

Hence, since 7^ 0, is a multiplicatively separated solution of (j2.9j) if and only 
if u is a separated solutions of (j2.5j) . □ 
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Remark 2. Although Eqs. ()2.8|) and ()2.9j) admit the same solutions for any 
choice of R 7^ 0, they are deeply different from the point of view of separabil- 
ity: equation (|2.9J) is treated in a natural way as an eigenvalue problem and its 
separation for all the eigenvalues j^E leads to /^-separation theory as developed 
by Moon and Spencer ^Hl QHj and by Kalnins and Miller [HI El Ell HTj On the 
contrary, in ()2.8j) the energy constant E is multiplied by the function R and E 
cannot be considered as a differential operator eigenvalue: in (|2.8jl E is a fixed 
parameter incorporated into the differential operator and the solutions are the 
eigenf unctions of the null eigenvalue only. Its separation must be treated in a 
way similar to the separation for Laplace equation ^3 Ej • The difference between 
these two view points is subtle but relevant (see also [3] about the correspondence 
between the quantum-mechanical description of hydrogen atom and the harmonic 
oscillators). Of course, also Eq. ()2.9|) can be read as a null eigenvalue problem 
by including the constant E into the potential function Ue = jkE — U. We 
adopt this last interpretation of (J2.9|) to define Fixed Energy R-separation (FER- 
separation) in Section 0] In this approach the function R plays the same role of 
the potential V and the metric tensor (g 1 -'): Eq. (|2.9jl is a PDE defined by G, V, 
R and E. By imposing that ()2.9|) admits a complete family of separated solutions 
0, in a meaning we will make precise later, we determine differential conditions 
on G, V } R and E. 

Remark 3. The change of unknown function (j2.1j) relates multiplicative and 
additive separation of any PDE, both in the case in which we are dealing with a 
single null PDE or with an eigenvalue problem. In the resulting equation (|2.5j) the 
parameter E admits two interpretations as in equation (|2.9|) : it can be considered 
as an internal parameter of the equation (FER-separation) or as an integration 
constant (classical R-separation). 

3 Geometry of the additive separation for a null PDE of 
order / 

We provide the geometric framework and the geometric characterization of the 
separation of variables (SOV) for a null PDE of arbitrary order /, extending the 
results of [4 (concerning the Hamilton- Jacobi case only) and of j2] (dealing with 
ordinary separation only). The main tool is the theory of separable connections 
as introduced by Benenti pQ . As in the case of the Hamilton- Jacobi equation with 
a fixed value of the energy, two different (but in a sense equivalent) definitions of 
separable solutions are given. Each of them is useful to prove different criteria of 
existence of a separable solution in a given coordinate system. 
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3.1 Separable connections 

Let M = QxZ^Qbe the trivial vectorial bundle with fiber Z over Q, where 
Q is a n-dimensional manifold with coordinates (q 1 ) (i = 1, . . . , n) and Z is a 
iV-dimensional vector space (over R or C) with coordinates (z A ) (A = 1, . . . , N). 

We call connection on M a regular distribution (in the Frobenius sense) C on 
M of rank n and transversal to the fibers: C : x C x C T X M with dimC^ = n 
and c x nr,z = {0}. 

A vector field X e X(M) is horizontal if X(M) C C(M) i.e., if X(x) G 
Locally a connection C is spanned by the n horizontal vector fields, called the 
generators of C 

d - 9 +r 9 

where CiA are functions on M called coefficients of the connection. A connection 
is integrable (in the Frobenius sense) if for any point P e M there exists a unique 
integral manifold (i.e., a n-dimensional immersed submanifold tangent to each 
Di) passing through P. By the Frobenius theorem and the particular form of the 
generators, we get that 

Proposition 4. A connection generated by the vector fields (Di) is integrable (in 
the Frobenius sense) if and only if [Di, Dj] = 0. 

If C is integrable, then there exists a foliation C of M made of integrable 
n-dimensional manifolds of C. Since C is a transversal distribution, the integral 
manifolds are locally described by za = Fa{o 1 , cb) where (cb) are real parameters 
(B — 1, . . . , iV)satisfying the completeness condition 



det 



dz/ 



7^0. 



dci 

The functions Fa are solutions of the first-order differential system of PDE on Q 

dz A 



dq i 

It follows that 



= C iA (q\ za). 



Proposition 5. The connection generated by the Di is integrable if and only if 
the corresponding system of PDE is completely integrable. 

Let S be a iV — 1-dimensional submanifold of M. We say that a distribution 
C is reducible on S if when it is restricted to the points of S it gives rise to a 
distribution Co on S. We have that 

Proposition 6. A distribution is reducible to a submanifold S if and only if its 
generators Di are tangent to S. 
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If F is a function on M such that equation F = implicitly defines the sub- 
manifold S, then C is reducible to F = if and only if DiF\ s = or, equivalently, 
DiF = AF for a suitable function A on M. 

Remark 7. If a connection C is integrable and reducible as a distribution to 
S, then the foliation £ = (L CA ) on M made of integral manifolds of C can be 
restricted to a foliation £$ on S made of integral manifolds of C, whose labels 
(ca) belong to a suitable iV — 1-dimensional submanifold V of R N . Up to a 
new parametrization of the (ca), we can assume without loss of generality that 
V is defined by = and £5 is described by the N — 1 parameters (c a ) = 

(Ci, . . . , Cjv-l). 



3.2 Equivalent definitions of additive SOV for a null PDE 

Inspired by the approach to the additive separation of Kalnins and Miller exposed 
in POj for PDE of the kind 7i = const., we extend the theory in a geometrical 
way to a single PDE of order I on Q 



H(q l , ...,q n ,u, u h Uij, uy.„&) = 



(3.1) 



where U; 



dfjii . . . are the derivatives (up to order I) of the un- 



known function u = u(q l ), with respect to the coordinates, with d{ = d/dq l . An 
(additive) separated solution of (jH.lj) is a function of the form 



u 



i=i 



where Si depends on the single variable q % satisfying (j3.1j) . If we look for additive 
separated solutions only, Eq. (|3.1|) reduces to 

H{q X , . . . ,q n ,u,u h uf \ . . . ,uf) 



0. 



where uf^ = uu and = for i ^ j. 

Definition 8. Given a PDE ()3.1|) of order I, the function -u 7 is a complete internal 
separated solution (internal solution) if: i) u 1 is of the form 



u 1 = ^Siiq'^Ca), 



a 



nl; 



8=1 



ii) for all (c a ) in a suitable open set of M. nl u 1 is a separated solution of (j3.1|) i.e. 



q n ,u,u I i ,u I ii ,. 



■, u u...ij 



iii) the following completeness condition holds 



rank 



du 1 

dc a 



dc n 



dc n 



du 1 . 



dc n 



1 - 



nl. 
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Definition 9. Given the PDE ()3.1|) of order /, the function u E is a complete 
extended separated solution (extended solution) if: i) u E is of the form 



u 



1, . . . ,nl + 1, 



ii) u E satisfies (j3.1j) . that is T~t(q l , . . . ,q n , u E , u E , u E , . . . , u E A ) = 0, for any (c a ) 
belonging to a suitable n/-dimensional submanifold of 



or, up to a trans- 
formation of the (c a ), for c n ; + i = 0; iii) u E satisfies the following completeness 
condition for all admissible values of q % and c a : 



det 



du 1 
dc n 



duf 
dc n 



9uj 

3Ca 



duj 

dCn 



^0. 



Remark 10. An internal complete solution defines a foliation of the submanifold 
7i = in n-dimensional leaves transversal to the fibers of the bundle M, via 
equations 



u = u 



Each submanifold is parametrized by the value of the nl parameters (c a ). Con- 
versely, an extended solution defines a foliation of an open subset of M in n- 
dimensional leaves via equations 



u = u 



dm 



U: 



(2) 



u 



(I) 



d ;U E . 



Each submanifold is parametrized by the value of the nl parameters (c a ). The 
foliation is compatible with the submanifold Ti, = in the sense that it is reducible 
to a foliation of 7i = 0. 

Proposition 11. In the coordinates (q 1 ) on Q, there exists an internal solution 
of the PDE ij) if and only if there exists an extended solution of (jJ. i|) in the 

(q l ) ■ 

Proof. Let u E [(f \ c a ) be an extended solution of (j3.1|) . Then, up to a reparametri- 
sation of the (c a ), we can assume that for all (c a ) with c n ; + i = the function 
u E is a solution of (|3.1|) . Hence, M 7 (g*,c a ) = u E (q\ Ci, . . . , c n /, 0) is an internal 
solution for (j3.1|) in the same coordinates. Conversely, let u 1 ^ 1 , c a ) be an internal 
solution of (j3.1|) . Let C m be the first column in the matrix 



du 1 
dc n 



dc n 



dc n 



du 



i 



OCn 



such that the n/-th order minor obtained by eliminating this column is differ- 
ent from zero. Such an index m exists because of the completeness condition. We 
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consider the two cases m = 1 and m > 1 separately. If m 
is a solution of (|3.1|) for all (c a ) such that c n ; + i 



1, then u E = u 1 +c n i + \ 



0; moreover, we have that 



M 



du E 
dc a 



dul 
dc a 



dc a 



du 1 



H 

dCa 













and the completeness condition holds. If m > 1, there exist unique q, r e Z such 
that m — 2 = qn + r with < r < n. Thus, the eliminated column contains the 
partial derivatives w.r.t. (c a ) of ra^fps with = g + 1 and j = r + 1. Then, 



it 



W + C ni+1 



is a solution of (|3.1|) for all (c a ) such that c ni+1 = 0. Moreover, the rank of 



M 



du 1 
dc n 



dc n 



dc„ 



dc n 



is maximal. Indeed, in the developing of det M E along the last row, the first m — 1 
addenda are null since M[ = for h < m, as well as those from the (m+ l)-th to 
the last, because the elements of the last row of M E with column index greater 
than m are all null. Hence, being det M E = (M E )^ +1 M^ = 1 • ^ 0, the 
completeness condition holds. □ 

Definition 12. We say that the null equation 7i = is separable in the coordi- 
nates (q l ) if it admits an extended solution or, equivalently, an internal one. 

Remark 13. In comparison with the geometric framework introduced in [3] 
to deal with a null Hamilton- Jacobi equation, the present case presents several 
differences. First of all, it is not possible to define complete internal or extended 
solutions without assuming that they are separated; moreover there is not a 
unique formula transforming an internal into an extended solution. This is a 
consequence of the fact that the coordinates in the fiber are not all of the same 
kind, but correspond to derivatives of u of different orders. 



3.3 Criteria for additive separation of a null equation 

Let us consider fibered coordinates (q\ za) on M = Q x Z — > Q, where (za) = 
(u, Ui, uf \ . . . , uf^), with i = 1, . . . , n, A = 1, . . . , nl + 1. Let 7i be a smooth 
function on M such that 7i = is a ^/-dimensional submanifold of M. We apply 
to this case the equivalence of integrability of distributions on vector bundles 
and complete integrability of normal first-order differential systems, recalled in 
section l3~Tl The following propositions provide necessary and sufficient conditions 
for the existence of an extended and an internal solution, respectively (see also 
J3] and references therein). 
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Proposition 14. In given coordinates (q l ) on Q, the following statements are 
equivalent: (i) there exists an extended solution u E of TC = 0/ (ii) there exist n 
functions Ri on M such that the distribution A generated by 

A = d t + Ui d u + ufd/dui + ... + u®d/du?- 1) + R l d/duf ) (3.2) 

is integrable and reducible on the submanifold TC = 0, that is 



DiRj = 0, V j ^ % 
DiH\n=o = 0. 



(3.3) 



Proof, (i) =>- (ii): Let us assume that a solution u E = £\ Si(q\ Ci, . . . , c/v) of 
TC — exists. By the completeness condition, (q\ci, . . . ,c^) are (non fibered) 
coordinates on M and the n functions on M 

d l+l 

Ri = WW s ^ ql,Cl '"' C7v) 

are well-defined. Going back the original coordinates we get functions 

Ri = Ri(q\ u, u i} uf\ . . . u\ l) ), 

and by them we construct the vector fields Di (J3.2j) . The distribution A generated 
by the fields Di is integrable because the corresponding first order system on M 



diU = Ui 
diUj = 5ijU 



(2) 



(3.4) 



diuf = SijRj 



is completely integrable. Indeed, 

u = u E = T,iSi{q\c A ), Ui = ±Si(q\c A ) ••• uf } = j^yS^, c A ) 

is a complete solution of the PDEs system, which describes a foliation C( CA ) = 
(L CA ) made of integral manifolds of A. Moreover, since u E is a solution of TC = 
for any (c A ) with c n i + i = 0, we have TC(q\ u E , . . .) = 0, that is every point of L CA 
belongs to TC = 0. Thus, the generators Di (by construction tangent to every leaf 
of C CA ) are tangent to TC = 0. Hence, A is an integrable distribution reducible 
to the submanifold TC = i.e., conditions ()3.3|) hold, (ii) =^> (i). We assume 
that there exist n functions Ri on M such that the distribution A generated by 
(|3.2|) is integrable and reducible on the submanifold TC = 0, that is (|3.2|) satisfy 
(|3.3|) . Hence, the first order PDE system associated with A (|3.4j) is completely 
integrable, that is there exist nl + 1 functions f A (q\ cb) such that z A = f A (q\ cb) 
are a complete solution of ()3.4|) that is satisfy ()3.4|) and the matrix (|^-) has 
maximal rank. The form of ()3.4|) implies that the solution is of separated type: 

u = J2 l S l (q\c B ), Ul = ^Si(q\c B ), ... uf ] = J £ 7 S i (q\c B ) 
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Thus, u E = J2 i Si(q l , cb) is an extended solution of TL = 0. Indeed, by con- 
struction it satisfies the completeness condition. Moreover, being the generators 
Di ()3.2j) tangent to Ti = by assumption, A is reducible to a distribution on 
TC — 0. Then, for any P in Ti = the connected leaf £( CA ) passing through 
P is entirely contained in Ti = 0. Let us consider the submanifold of R ni+1 
V = {(eg) I L [cb) cH = 0}. Then, for all (c B ) G V u E satisfies H = 0. □ 

Proposition 15. In a gwen coordinate system (q 1 ) on Q the following statements 
are equivalent: (i) there exists an internal solution u 1 = £V Si(q\ c a ) of7i = 0; 
(ii) there exist n functions Ri on M such that the distribution A generated by 
the vector fields flJ.^j) is reducible to a distribution A on the submanifold Ti = 
and Ao is integrable, that is 

DiH\n=o = 0, ,„ j.\ 

A^k=o = o, Vj^i. 1 ' 

Proof, (i) =>- (ii): if there exists an internal solution, then there is also an ex- 
tended solution in the same coordinates. Hence, by Proposition El there are 
functions Ri satisfying ()3.3|) and in particular ()3.5|) . (ii) =>- (i): Being A in- 
tegrable, the manifold Ti = is foliated by integral manifolds of A locally 
described by solutions of the PDE system (J3.4|) depending on nl parameters (c a ) 
with separated form 

U = Tli Si(q l , C a ), Ui = ■£jiSi(q' 1 , c a ), ... u\ — r^y Si(q l , c a ) 

satisfying Ti = 0. Thus, u 1 = ^ Si(q\ c a ) is an internal solution of Ti = 0. □ 

Theorem 16. The separation of a null equation Ti = occurs in a given coordi- 
nate system if and only if one of the following equivalent conditions holds: 

1) there exist n functions Xi(q, u, Ui, .., ') such that the operators (jff.jjj) satisfy 

D{H = XiTi., , „s 

DiRj = 0, Vj ^ i. {6 - 0> 

2) there exists (locally) a function A(q l , u, Ui, uf ) such that the equation 
Ti/A = h is separable in the ordinary sense, i.e. the operators satisfy 

A(W/A) = 0, 

3) the operators (j, 1 ?.^) satisfy 



DiH = 0, , , 

A^|«=o = 0, Vj^i. ^ 
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Proof. Condition 1) follows straightforward from (|3.3|) by applying Hadamard's 
Lemma: F\ n=Q = <^ F = XH. Condition 2) is equivalent to condition 1). 
Indeed, from (|3.3jl i we get A, = A logTi. Since A-A = DjD iy we have AAj = 
DjXi and locally there exists a function A such that 

Xi = Di log A. 

Thus, condition ()3.3)l i becomes DihiTi = Di In A and, being Di linear, 

A(W/A) = 0. 

Hence, condition 1) holds iff there exists A satisfying (|3.7jl . Condition 3) is a 
consequence of Proposition HH1 Indeed, by applying Hadamard's Lemma to ()3.5j) 
and (|3.8|h we get 

3 W' WI { Ai£=^«, V^i, (3 - 9) 



3 (i/fefe) 



(3.10) 



AH = o, 

respectively. We show that ()3.9|) is equivalent to ()3.10j) . Clearly, (|3.10|) implies 
(|3.9J) by choosing Aj = 0. Conversely, we assume that (|3.9J) holds and we denote 
by the operators (|3.2jl and the associated functions i?j satisfying AH = 0. 

moreover we denote by Df, Rf those satisfying D{H = XTi. Then, we have 



D9 = A A - TCh. 



d 



du 



(0 ' 



IX; — IX; — J^ilX. 



where Ti = (-^jjj) and by a straightforward calculation we get 



DfR^ = a u #, u + U 



D 



Thus, since D^Rj = fiijTC, there exist functions i/y such that D^R® = i^W. □ 

Remark 17. The n additional unknown functions (Aj) play the role of La- 
grangian multipliers. It is remarkable the fact that the n functions (Aj) can 
be replaced by a single unknown function A. 

Remark 18. We recall that the conditions for the separation in the ordinary 
sense of the equation TC = h are lU HHj 



D{H -- 
DiRj 



0. 

= o, 



(Vj + i). 



(3.11) 
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4 Definition and criteria for F£"i?-separation 



By Proposition ^ we can extend the notion of completeness for separated so- 
lutions of (|2.5jl to FER-separable solutions of (jl.ljl in a natural way. However, 
since Ti does not depend on u, but only on its derivatives, the solution of (J2.5j) 
is defined up to an inessential additive constant (corresponding to a constant 
factor for ijj). Hence, in this case we can consider without loss of generality only 
the 2n functions (za) = uu) as unknowns (see Remark 2.4 in |2j). It follows 
that in the generators Di the term Uid u can be disregarded. Moreover, being the 
unknown functions 2n, the number of parameters entering in a complete solution 
of the null PDE ()2.5|) is 2n — 1, according to Section f3.ll Therefore, we define 
FER-separation as: 

Definition 19. The Schrodinger equation (jl.lj) is FER-separable in the coordi- 
nates (q % ) - called FER-separable coordinates for (jl.ljl - if there exist a value of 
the energy E and a function R such that Eq. (jl.lj) . admits a solution of the form 



^ = R\\(j)i{q\c a ) 



(4.1) 



where c a are 2n 



1 parameters satisfying the completeness condition 



rank 



dc a 

dvj 
dc n . 



2n- 1, 



(4.2) 



Remark 20. The meaning of completeness is the following: given R and E 
allowing FER-separation, for any choice of 2n — 1 numbers (c a ) = (pi,kj), i = 
l...n, j = 1, ...,n— 1, and of a point go £ Q, then there exists a unique 
separated solution of (J2.9j) such that (^ ^ L ) 90 = (/ii, fcj). Moreover, we observe 
that in Definition there are no assumptions on the form of the functions (pi 
(as for instance fa = e biql with bi G M) as in the case of free separation of the 
Schrodinger equation [2j. 

We apply the third condition of Theorem to study the separation of Eq. 
(12. 5 j) that is the null equation Ji = where 

H = g ij u lUj + g u u ti - Fin + ^E-U. (4.3) 

We assume here and in the following that g u ^ for all i. This assumption is 
obviously fulfilled for a proper Riemannian manifold. This means that we are 
excluding null coordinates. In our case, we have I = 2 and Ti independent of u. 
Hence, both the operators (J3.2j) and the associated functions Ri do not depend 
on u: 

D i = d i + u ii — + R i —. (4.4) 
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Since D/H = 0, the functions Ri have the form 



Ri(q\ui,Ui 



d{H + u u —— 

OUi 



du 

We want to find conditions for which DiRj\w = Q = (i ^ j 



(4.5) 



Lemma 21. The function DiRj is a polynomial of second degree in the variables 
(tin, . . . , u nn ), whose quadratic part is 

gij 

--— "»»xr 1,3 n.s. i^j. (4.6) 
9 

Proof. Recalling that J^- = g n , by inserting (j4.3|) in ()4.5|) . we get 

Rj = i[ -2^ h u i5 u h - d 3 g hh u hh - d,g hk u h u k + d,T h u h + Pu n + 8 3 U}. 

Since i?j and all its partial derivatives contain terms at most linear in the variables 
(u hh ) and being 



dR dg" 1 

DjRj = dsRi + Uu— — 3 . . Rj 

3 3 dui g3J 



a straightforward calculation shows that the only term of DiRj which is quadratic 
in (u hh ) is (USD. □ 

Proposition 22. The coordinates (q l ) allowing FER- separation i.e, separation 
of the null equation 7i = (j£.<5j) for a suitable R, are necessarily orthogonal. 

Proof. By condition 3) of Theorem EE the null equation 7i = is separable if 
and only if for all i ^ j, the functions DiRj vanish on the surface TL — 0. Since 
the function ()4.3j) is linear (uhh), the equation TC = is equivalent to 



9 



4? ( + 9 hk u h u k - f V + ^E-U). 



Hence, inserting the expression of Uu into DiRj, we get (by Lemma ETTj) a quadratic 
polynomial in the (u aa ) (a ^ i) with quadratic part 



9 13 9 aa 

a* 9 9 



which vanishes if and only if g ij = for i ^ j. □ 

Remark 23. By Propositionl221coordinates allowing FER-separation of Eq. (jl.lj) 

are necessarily orthogonal. However, nonorthogonal .R-separation (see [HIE]) is 
possible by imposing constraints on the form of some factors (f)f this case, which 
has been described in j2] for ordinary separation of Schrodinger equation and 
called reduced separation, is not examined in the present paper. 
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In orthogonal coordinates, (|4.3|) becomes 

L W 



H = g ll (u tt + v%) - Vu t + ±E-U (4.7) 



and the functions ()4.5|) 

R, = ^[-2g^u j3 ^ - d j9 hh (u hh + ul) + djt h u h + Pu n + djV). (4.8) 
9 

In order to write in a shorter way DiRj we introduce the following notion 

Definition 24. A Stackel operator is the second order differential operator 
such that for any f:Q—*M. 

SM = dy - d j \n(g ,; )i) : f - ih H9 jj )9jf. (4.9) 

Proposition 25. Equation is separable in the orthogonal coordinates (q l ) 

if and only if 

djV - rdj In g" = 0, (4.10) 

Sij{9 hh ) S tJ (g kk " 



0, Vh,k, (4.11) 



S ij (U)g hh -S ij (g hh )(U-£E) = 0, V /i. (4.12) 
Proof. By applying operators ()4.4|) to the functions (|4.8|) and using (|4.9j) we get 

A-Rj = 4?[(^ f * - '"^ 111 //" )"" + - f^i ln^>^ 

-^•(/ ft )K A + t#) + Sy(f> h + S^(C/)]. (4.13) 

We impose that the function -Dj-R., vanishes on the surface 7i = that is for all 
values of (uh,ul,u aa ) with a^i and for 



Uu = uf 



^ ( E ^ (Uaa + <) - t h U h + &E-U). 



9 a+i 



By inserting the expression for ua in ()4.13|) we see that the only coefficient of u\ 
is djT 1 — T % dj \ag 11 . Thus, ()4.1U|) must hold for any couple of indices i ^ j. Then, 
assuming (|4.10|) . DiRj\-H=o becomes (up to a factor gjj) 

Sij(g")„ a g o f„aa\\„ j_fc (rh\ M£!2r*\. S i j(g")(U--2 I E) 



y-ir - SiA9 aa )\ v a + U,-(r fc ) - ^-v h j u h - ^ ; ;„ ^ + s^u) = o 

with v a = u aa + v? a . By equating to zero the coefficients of v a , Uh, and of the 
0-th order term, we get two conditions equivalent to (J4.ll)) . (|4.12j) and moreover 

S ij (r h )-**p-t h = 0. 

However, this last condition is disregarded being a consequence of (|4.10ff4~TT|) . 
Indeed, by (fOB we get S^T* 1 ) = f 1 ^^)^. □ 
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Remark 26. If U - Jr£ = 0, i.e., if R satisfies (JET]), then fEEj) is always 
satisfied. 



Proposition 27. If (q l ) are separable coordinates for \2. 5|) . i/ien the function R 
is locally determinated by integrating equations 

2d i \nR = T l -Uq 1 )- (4.14) 

Proof. For f * ^ 0, we have that (|4.1(J|) is equivalent to dj (P /g n ) = for all j ^ i 
i.e., to P = g ll ^i{q l ) {i n.s.). By definition of modified contracted Christoffel 
symbol (|2.7jl . since the coordinates are orthogonal, we get Tj — 2<9j lni? = £i(g l ), 
that is ()4.14j) . The differential system is integrable in orthogonal coordinates 
satisfying (|4.11|) . Indeed, a straightforward calculation shows that 

djTi = O.Y, <> aj In ^ = d% In / •<:> = 

□ 

Remark 28. By using ()4.14|) we can write AR/R (included in the modified 
potential U) in terms of the contracted Christoffel symbols Tj as 

AD 1 

— = i /(2^r, - r, 2 + c 2 - 2d&). (4.15) 

We shortly recall some definitions arising from the theory of separation of 
geodesic Hamilton- Jacob i equation (see |lj). 

Definition 29. Orthogonal coordinates (q l ) are separable if the metric g n is a 
Stdckel metric i.e. it is a row of the inverse of a regular matrix S = (y^ j) ) (Stdckel 
matrix) whose elements depend on the coordinate corresponding to the lower 
index only. Orthogonal coordinates (q l ) are conformal separable if the metric g u 
is a conformal Stdckel metric i.e., there exists a function A such that g u /A is a 
Stackel metric. A function / is a pseudo- Stdckel factor (resp. Stdckel factor), if 
it can be written as f — g u <pi{q l ), where (g tl ) is a conformal Stackel metric (resp. 
Stackel metric). 

Theorem 30. Necessary and sufficient conditions for the FER-separation of 

in a given coordinate system (q l ) (under the assumption g u ^ are: 
(1) the coordinates are orthogonal; (2) the coordinates are conformal separable; 
(3) the function 

^{E-V) + \g ii {2d i V i -VD 
is a pseudo-Stdckel factor. In this case R is any solution of equations OTT^J). 



Proof. By Proposition |2U we have that coordinates are necessarily orthogonal. 
Moreover, by Proposition |2H1 i n orthogonal coordinates FER-separation occurs 
if and only if conditions ()4.10|) . (j4.11j) and ()4.12j) hold. Equations (|4.11j) means 
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that the orthogonal coordinates q % are conformal separable coordinates, condition 
()4.12j) means that U—j^E is a pseudo-Stackel factor (see jlj) and, by Proposition 
E3 Eq. (j4.14|) determines R up to separated factors. By (|4.15J) . we have that we 
can replace U — j^E by the function m{E — V) + \g ll (diFi — ^Tf), disregarding 
the term \g n 9i which is a pseudo-Stackel factor for any choice of £j in (|4.14|) . □ 

Remark 31. In the ordinary separation for Schrodinger equation, condition (3) 
is split into the compatibility condition for the potential (V = g u 0i(q 1 )) and 
the Robertson condition (diTj = 0). Moreover, the coordinates are in this case 
necessarily separable, while in the case of FER-separation conformal separable 
coordinates are also allowed. 

Remark 32. For the Laplace equation, condition (3) of Theorem I3TH means that 
\g hh {2dhTh ~ T^) is a pseudo-Stackel factor, that is 

„hh 

I y —(2d h T h - I*) = g%tf). (4.16) 

We remark that this condition in general does not imply (see [B]), 

2 3^,-1* = 40,(0*). Vz (4.17) 

For example toroidal coordinates (see Section EJ) satisfy ()4.16j) since g hh {2dfXh — 
T\) = g 33 , but they do not satisfy (|4.17jl . Moreover, even if the (q l ) are separable, 
()4.16|) does not imply the Robertson condition diTj = (by (|4.14jl Robertson 
condition implies that R is a separated function and /^-separation is then called 
trivial JD])- A counterexample of that is given by the conformally flat metric 
allowing non trivial .R-separation in jS]. In ^3] condition (J4.17|) is linked to the 
fact that symmetry operators can be put in the reduced form. 

Proposition 33. If (q l ) are FER- separable coordinates for (jl.l)) for two values 
of the energy E\ ^ E2, then (q l ) are separable and allow R-separation for all 
values of the energy. 

Proof. If (jj~I2D holds for both E x and E 2 , then Sij(g kk ) = 0, k = 1 . . . n and 
QUE)) hold for all Eel. □ 



4.1 Separated equations 

If q l are conformal separable coordinates for the orthogonal metric G = (g u ), then 
for any pseudo-Stackel factor / we have that (g u /f) is a Stackel metric. Hence, in 
coordinates allowing FER-separation we have that the conformal metric G with 
components 

9" = -ry^-r (4-18) 
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is a Stackel metric and that Eq. ()2.5|) is equivalent to 



r(««+«i-&«t) = -i, (4-i9) 

for any choice of the functions £,i{q 1 )- In order to integrate Eq. (j4.19j) bv separation 
of variables, it should be considered as the n-th equation of the system 

+ u\ - itUi) = cij, (a n = -1) 

where (pfy is the inverse of a Stackel matrix associated with the Stackel metric G 
and (flj, . . . , a n ) are real constants. By inverting this system, we get the separated 
equations 

n-l 

uu + u] - £ iUi = «^1 Q) - ^ , (4-20) 

a=l 

corresponding by using (|2.1()j) to the Riccati equations for the separated factors 

<f>i 

$ -M- <P? a, & = 0, (a n = -1). (4.21) 

Being R locally determined (up to separated factors) by Eq. (|4.14jh it is 
evident that we always can choose £j = for i — 1, . . . , n so that 1)4.14)1 becomes 

2dMR = Ti. (4.22) 

It follows that 

Proposition 34. In FER-separable coordinates, by choosing R satisfying ( . 20$ , 
separated equations ( 4-2(Jty , \4-%l\ have the canonical form 

2 (i) 

Remark 35. Even if they do not explicitly appear in the separated equations, 
the potential V and the value E are contained in the components Stackel metric 

Remark 36. The above described method of separating the variables gives rise 
to two kinds of constants in the separated solutions: the n — 1 constants dj (j = 
1, . . . , n — 1) of Eq. ()4.21jl that we call separation constants and the integration 
constants of Eq. (|4.21|) . Even if the 2nd order Eq. (|4.21|) give rise to 2n constants, 
n of them are inessential constant factors of the 0j, irrelevant for our purposes as 
explained at the beginning of the Section. Therefore, there are only n non-trivial 
integration constants {pi = <j>i) qo ) which, together with the n — 1 (cij), form the 
set of 2n — 1 constants satisfying the completeness condition fj4.2jl . The relation 
between the constants (a,-, bj) and the constants (foj, kj) introduced in Remark I2TH 
with [kj) = {</>"/ <j>i) qo ) is given by (|4.21|) evalueted in the initial point q . 
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In the following subsections, we shall see that, in analogy with the ordinary 
separation of the Schrodinger equation (see [TUl 13 E| and references therein), 
the separation constants are in fact the eigenvalues of n second-order pairwise- 
commuting conformal symmetry operators of the equation related to quadratic 
first-integrals in involution of an associated Hamilton- Jacobi equation. 

4.2 Conformal Killing Tensors and FER-separation 

Fixed-Energy separation of (|2.5jl is characterized by the existence of CKTs of 
the metric G in conformal involution with common eigenvectors (see 0J). In- 
deed, conformal separable coordinates are associated with n symmetric 2-tensors 
Ki, . . . K n = G which are Killing tensors with respect to the metric G fl4.18|) . 
simultaneously diagonalized in the (g*), and with contravariant components 

Kf = ifi\ j} , K} h = 0, forz^/i. (4.23) 

Let us consider the cotangent bundle T*Q with canonical coordinates (q l ,Pi) and 
Poisson brackets of functions A, B on T*Q defined by 

4-1 dpi dq l dq i dpi ' 
Let us construct the quadratic polynomials in the momenta associated with (Kj) 

P Kh = Klp iVi h^n, P Kn = gV PiPj . (4.24) 
We have that OH the functions P Kh defined by and flPSj) satisfy for h, 

j — 1 TL\ 

{P Kh ,P K] } = 0, (4.25) 

i.e., Pk k are pairwise in involution first integrals of the geodesic Hamiltonian Pk„, 
and the geodesic Hamilton- Jacobi equation g^PiPj = h is separable in the (q z ). 
Hence, we get 

Proposition 37. An intrinsic necessary condition for the existence of FER- 
separable coordinates for the Schrodinger equation (ji.ip is that there exist n — 1 
CKTs Kx,...K n _x such that (i) Ki,...K n _i,G are linearly independent, (ii) 
Ki, . . . K n _i have common eigenvectors, and (Hi) Eq. Q-^.i?<5| ) is satisfied. 

4.3 Symmetry operators 

Starting from separated equations (J4.2(Jj) and (|4.21j) we build conformal symmetry 
operators of (jl.lj) . Interpreting the left-hand side of (J4.2(Jj) as the results of n 
differential operators acting on u, 

u ^ (p^un + uj - £iUi), (4.26) 
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and inserting ()2.10|) . (|4.14jh ()2.2|) in (|4.26|h we get n operators Hj on if) defined 
by 

Hjij) = Ajtp - ^AjRtp, (4.27) 
R 

where 

A.j : if) i— > - IW) = 4)^, 8i:il>»cfol>- mip. (4.28) 

For j = n we have A n = A/(-jkE — U) and for a n = —1, i7 n ^> = a n ?/> is equiva- 
lent to the Schrodinger equation (jl.l|) . By Proposition |23 since f)j\/g u are the 
eigenvalues A* w.r. to G of n diagonalized KTs of G, defined by 

and since the A^ satisfy the following intrinsic Killing-Eisenhart equations (see 

EHHi) 

^ = (A}-Aj)^log^, 
a direct computation shows that 

Proposition 38. Let (q l ) be orthogonal conformal separable coordinates and (K,) 
be n independent KT of the conformal metric G ft4-18\) in involution and simul- 
taneously diagonalized [K 1 - 1 = 8 %h (p^), with K„ = G. Then, the operators Hi 
pairwise commute i.e., for j, k — 1, . . . , n 

[Hj, H k ] = HjH k - H k Hj = 0. (4.29) 

By expanding [Hi,S], where S denotes the Schodinger operator S : ip i— > 
-Aif) + £V, by KT7\i and KM . we see that H^ are conformal symmetry oper- 
ators 9\ for S. It follows that 

Proposition 39. FER-separable solutions if) are common eigenvectors of com- 
muting second order conformal symmetry operators for the Schodinger operator 
S. 



5 An example: toroidal coordinates on the Eu- 
clidean three-space 

Let us consider toroidal coordinates (q l ) = (rj, 9, ip) on the Euclidean three-space 
E 3 . These coordinates have been applied also to biophysical systems jT3J. The 
transformations to Cartesian coordinates (x, y, z) are 

a sinh rj cos ip 

cosh 77— cos 6 ' 
= a sinh r, sin g R+ 

a cosh rj— cos 8 ' 

a sin 8 

cosh rj— cos 8 ' 
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The coordinates hypersurfaces q l = const, are toroids, spherical bowls, and half- 
planes through the z-axis, respectively |17j . The non vanishing contravariant 
components of the Euclidean metric in the (q l ) are 



n 22 (cosh r] — cos ( 



9 



33 



(cosh i] — cos ( 



a" a 2 sinh 2 1] 

We consider the FER-separation of Schrodinger equation for E = 

- ±ViP = 0. 



(5.1; 



Since for toroidal coordinates we have g hh {2dhX'h — T|) = g 33 , which is a Stackel 
factor, condition (3) of Theorem EDI provides the form of all possible potentials 
allowing FER-separation of (|5.1|) 



V = (cosh?? - cos6) 2 fM + f 2 (9) + -—^fs(<f) 

\ sinh r\ 

where are arbitrary functions. By integrating (j4.14j) we get 



/ cosh rj — cos 9 \ 2 
\ sinh r\ 



(5.2) 



for all & = (canonical separated equations), while we get R = (cosh rj — cos 6)? 
for the choice £i = cotanhr/ (non canonical separated equations)) The conformal 
metric g %% = g n jo with 

a = (cosh r? - costf) 2 (fi{rj) + f 2 {9) + (sinh 77) ~ 2 / 3 (</?)) 

is a Stackel metric and associated with the Stackel matrix (depending on V) 



S 



The separated equations are 
d 2 



fi — 1 — sinh 2 77 

h 1 

/a 1 



, ^ + cotanhr/— — + (/j - c 2 - c 3 (sinli77) 2 )0i = 
a?7 z drj 



d 2 < 



dd 2 



d 2 < 



d V 2 



+ (/a + c 2 )0 2 = 

+ (/3 + C 3 )03 = 
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where the term in -^p disappears if R is chosen to be (|5.2|) . Two conformal Killing 
tensors in involution associated with toroidal coordinates are 

K 1 = d e ®d e + f 2 G K 2 = d v ® <9 V + / 3 G 

whose components in Cartesian coordinates (disregarding the term proportional 
to the metric tensor) are 



1 

2a? 



2x 2 z 2 2xyz 2 
2z 2 y 2 



-(x 2 + y 2 + z 2 )xz 
-(x 2 + y 2 + -2 2 )?/2; 
(x 2 + y 2 + -z 2 ) 2 



-xy 




-xy 
x 2 




The research of non trivial examples of FER-separation for a single value E ^ 
in dimension n > 2 is in progress. 
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